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1. Complex Parameters in the MSSM 
Higgs Sector 

The Higgs sector of the Minimal Supersym- 
metric Standard Model with complex parameters 
(cMSSM) consists of two Higgs doublets 
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which form the following Higgs potential: 
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has the form 
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where mfi,H,A denote the tree-level Higgs masses 
and the self-energies' momentum dependence has 
only been suppressed for compact notation. It 
should be noted that in general M^ is symmetric 
but not Hermitian. In the approximation of van- 
ishing external momentum (g^ = 0), one can ob- 
tain the higher-order corrected mass eigenstates 
via a unitary transforination from the tree-level 
states: 



^hA /Un U12 C/13 
^2 = U21 U22 U23 

M \U3i U32 c/33 




(4) 



The Higgs potential contains two complex phases 
£,, arg(mf2)- The phase arg(mf2) can be rotated 
away [112] and, at tree level, ^ has to vanish in or- 
der to fulfill the minimum condition of the Higgs 
potential, so there is no CP- violation at tree level 
and the spectrum contains five states of definite 
CP-parity: /i, H, A, iJ±. 

CP-violating effects are induced by complex pa- 
rameters that enter via loop corrections: the tri- 
linear couplings At.b,T, the Higgsino mass param- 
eter /x, and the gaugino mass parameters Mi ^2,3- 
They yield 'ShA, '^ha 7^ and induce mixing be- 
tween h, H, and A [5]. The Higgs mass matrix 



2. Higgs-boson Self-Energy Corrections in 
FeynHiggs 

2.1. Higgs-boson Masses 

The following contributions to the mass ma- 
trix and the chargcd-Higgs-boson self-energy are 
taken into account: 
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® Leading 0{atas) cMSSM two-loop corr. g]. 

@ Leading 0{a'^) and subleading 0{abas, atab, 
al) two-loop corrections evaluated in the real 
MSSM (rMSSM) [SMTT] . For the treatment of 
phases see Sect. 12. 2] below. 

® Full one-loop evaluation (all phases, g^ depen- 
dence) |14| and leading non-minimal flavour- 
violating (NMFV) corrections [TU] . 



FeynHiggs [11112113114] performs a numerical 
search for the complex roots of det M^((7^) which 
are denoted A^^ . We decompose Ai^ = M^ — 
iAir, where M is the mass of the particle and 
r its width, and then define the loopcorrected 
masses according to Alf^ ^ M/i^ ^ M^^- 

The Higgs masses are thus determined as the 
real parts of the complex poles of the propagator. 
Complex contributions to the Higgs mass matrix 
(from ImS) are taken into account [14|15j . 

2.2. Two-loop Corrections in the cMSSM 

The full phase dependence is taken into account 
in the complete one-loop [14j and 0{atas) two- 
loop [4] contributions to the Higgs self-energies. 
We used the approximation of vanishing external 
momenta and vanishing electroweak gauge cou- 
plings in the evaluation of all two-loop diagrams. 

The tlCplxApprox flag controls the treatment 
of phases in the part of the two-loop corrections 
known only in the rMSSM so far. The following 
values are possible: 

• = all corrections: 0{atas,abCts,ct^,cttC(b,C(b) 
in the rMSSM, 

• 1 = only the cMSSM 0{atas) corrections, 

• 2 = the cMSSM 0{atas) corrections combined 
with the remaining corrections in the rMSSM, 
truncated in the phases, 

• 3-6 = the cMSSM 0{atas) corrections com- 
bined with the remaining corrections in the 
rMSSM, interpolated in the phases. 

• New in 2.7 is the choice of interpolation vari- 
ables: At,Ab,M3,^l (3), At,Xb,M3,n (4), 
Xt,Ab,Ms,^i (5), Xt,Xb,M3,ti (6). 

FeynHiggs thus not only has the most pre- 
cise evaluation of the Higgs masses in the 



cMSSM available to date (using the Feyn- 
mandiagrammatic approach), but can also de- 
liver an estimate of the uncertainties due to the 
rMSSM-only parts, by comparing the output for 
different values of tlCplxApprox. 



2.3. Improvements in the Af, Resummation 

Corrections to the bottom- Yukawa coupling are 
potentially large and important phenomenologi- 
cally, in particular for large values of tan/?. 

A resummation of the terms 0(a" tan" /3) and 
0(a" tan" /3) can be performed with the help of 
Aft P9] . 

Some care has to be taken when applying the 
resummed Af, corrections in FeynHiggs because 
the Af, terms inserted into the one-loop Higgs self- 
energies have to be consistent with the two-loop 
corrections selected. Whenever the corrections @ 
in Eq. (|2.ip are enabled (tlCplxApprox ^ 1), we 
also have to use the Af, of Ref. |6j. The deter- 
mination of Af, was improved in this direction in 
Version 2.7, avoiding incomplete higher-order cor- 
rections at 0{abas) and beyond. 

In all other instances (e.g. decays, production 
cross-sections) , the Af, of [17] is used in FeynHiggs 
2.7, computed self-consistently according to the 
procedure in jl8] . 



3. Decay Rates and Production Cross- 
Sections 

3.1. ft-i — J> fjfk at One-Loop Precision 

The Higgs decays to fermions, hi — )• fj fk are 
now available at full one- loop precision E^- The 
real gluon (photon) which cancels the IR pole is 
treated fully inclusively j25j . 

The routine FHCouplings has an additional 
flag fast to switch off the computation of the 
off-diagonal decays {j ^ k) for efficiency. These 
decays hardly contribute to the total decay rate 
and can safely be neglected. 

The phenomenologically important resummed 
Af, corrections are still taken into account in hi -^ 
bb, with the corresponding one-loop contribution 
subtracted to prevent double counting. 
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3.2. Improvements in the gg -^ h Produc- 
tion Cross-Section 

For the Standard Model estimate, we use the 
NNLL prediction of Ref. [H]. 

The MSSM production cross-section is esti- 
mated in the effective-coupling approximation 
[23j by multiplying pieces of the amplitude with 
the corresponding state-of-the-art amplitude cor- 
rection factors as follows: 
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The four sectors are clockwise labelled LL, LR, 
RR, RL with entries 
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The factors cf^°, cf '^^^j Cf,_r, Cb^i are taken from 
Ref. [21] and by construction reproduce the top 
and bottom loops at NLO including the interfer- 
ence term. Cf is taken from Ref. 22 . The LO 
cross-sections are parameterized with ■mb{'mb)- 

The crS'^''^L° ig jjigQ ^j^j^gjj fj.Qj^ Ylei. [H] and 
contains only SM top and bottom loops. 



4. Non-Minimal Flavour Violation 

In the NMFV MSSM, the sfermion flavours are 
allowed to mix with each other, i.e. the mixing is 
6x6 rather than 2x2: 
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and likewise for the sdown sector. 

The mixing matrices R^''^ diagonalize the mass 
matrices (M^)"'^ = (M^jpy)"'^ + A"'^, where 



L,q 



Mq ,+m;+ cos 2/3 (r| - Q,sl.)ml , 
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Technical remark: FeynHiggs 2.7 keeps the MFV 
mixing matrices U exactly 'on top' of the NMFV 
matrices R, such that the MFV entries are au- 
tomatically updated whenever the NMFV entries 
change and vice versa. 

The most immediately notable effect comes 
from the LR(RL) sector, as the AJ^ enter the cou- 
plings directly, e.g. 



A 



^ [ rrid^, Rtg+sRtg, {5cjg' /i + A'pg tan /?) - 
«■«' md^Rf:gRlg,^^{5gg.^Ji* + Alg, tan/3) 



This enters the Higgs masses through the A self- 
energy and can lead to sizable effects. 

The main constraints on the A"'"^ come from 
low-energy observables. Currently included in 
FeynHiggs are 6 — )► 57 and AM^, both at one- 
loop level including NMFV effects, with more to 
follow. 

5. Output of FeynHiggs 2.7 

We give a short overview of the output routines 
of the FeynHiggs library. 

FHHiggsCorr - All Higgs-boson masses and mix- 
ings: Mhi,2,3^ Mh±, acff, UHiggs, ZHiggs. 

FHUncertainties - Uncertainties of the masses 
and mixings. 
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FHCouplings - Couplings and Branching Ratios 
for the following decay channels: 

hi,2^s^hfj,ll,ZZ^*\WW^*\gg, 
h^Z*,hihj,H+H-, 



r*>f,h,w^*Jf',x°xf, 



JJ± _^ f(*) fl 
t^W+b,H+b, 



plus the corresponding ones of an SM Higgs with 
mass Mh,: /i^Mg ^ /,/„ 77, ZZ(*\ WW^*\ gg. 

FHHiggsProd - Higgs production-channel cross- 
sections (SM total cross-sections multiplied with 
MSSM effective couplings, see Ref. |l9l23l27) ) 

• gg ^ hi ~ gluon fusion. 

• WW -^ hi, ZZ -^ hi - gauge-boson fusion. 

• W -^ Whi, Z -> Zhi - Higgs-strahlung. 

• bb ^- bbhi - bottom Yukawa process. 

• bb ^ bbhi - ditto, one b tagged. 

• tt ^ tihi - top Yukawa process. 

• tt ^ tthi - stop Yukawa process. 

FHConstraints - Electroweak precision observ- 
ables, see e.g. Ref. [26] for details: 

• Ap at 0{a,aas), including NMFV effects. 

• Mvy, sin^ Ocfi via SM formula + Ap. 

• BR(6 -> S7) and AM^ at one-loop level includ- 
ing NMFV effects [H]- 

• (5^1 — 2)suSY including full one- and leading/ 
subleading two-loop SUSY corrections. 

• EDMs of electron (Th) , neutron, Hg. 

6. Download and Build 

• Get the tar file from feynhiggs.de, unpack and 
configure: 

tar xfz FeynHiggs-2 .7 .0. tar .gz 
cd FeynHiggs-2 .7 .0 
. /configure 



• "make" builds the Fortran/C-t--|- part only, 
"make all" builds also the Mathematica part. 

• "make install" installs the package. 

• "make clean" removes unnecessary files. 

7. Usage 

FeynHiggs has four modes of operation: Li- 
brary mode. Command-line mode, Web mode, 
and Mathematica mode. 

7.1. Library Mode 

The FeynHiggs library libFH. a is a static For- 
tran 77 library. Its global symbols are prefixed 
with a unique identifier to minimize symbol colli- 
sions. The library contains only subroutines (no 
functions), so that no include files are needed (ex- 
cept for the couplings) and the invocation from 
C/C-l— I- is hassle-free. Detailed debugging out- 
put can be turned on at run time. All routines 
are described in detail in the API guide and on 
man-pages. 

7.2. Command-Line Mode 

The user submits a text file, such as 



MT 


172.6 


MSusy 


500 


MAO 


200 


Abs(M_2) 


200 


Abs(MUE) 


1000 


TB 


5 


Abs(Xt) 


1000 


Abs(M_3) 


800 



to the FeynHiggs executable with a command like 

FeynHiggs file Lflagsl 

where the flags are optional. The output is a 
human-readable version of the results. The table 
utility converts the output to machine-readable 
format, for example 

FeynHiggs file Iflags'i I table TB MhO > out 

Loops over parameter values (parameter scans) 
are possible: 

• MAO 200 350 50, linear: 200, 250, 300, 350, 

• TB 5 40 *2, logarithmic: 5, 10, 20, 40, 

• TB 5 50 /4, number of steps: 5, 20, 35, 50. 
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7.2.1. Command-Line Mode Scripted 

If a "-" is specified as file name in command- 
line mode, FeynHiggs reads the input from stdin. 
This allows to script FeynHiggs sessions, such 
that parameters, flags, and possibly environment 
variables are preserved in a compact shell script. 
For example: 



#! /bin/sh 








make 1 1 exit 1 






FHDEBUG=2 


/build/FeynHiggs 


- ${1: 


-400302103} « END 


MT 


173.1 






MSusy 


3000 






MAO 


1000 






Abs(M_2) 


2500 






Abs(MUE) 


2000 






TB 


5 






Abs(Xt) 


1000 






Abs(M_3) 


2000 






END 









make I I exit 1 updates the FeynHiggs exe- 
cutable and terminates with an error if any prob- 
lem occurred in the build. 

The prefixed FHDEBUG=2 sets an environment 
variable, here to increase the debugging level. 

There are default flags, overridable through 
command-line argument 1 (${1 : -400302103}). 

The actual parameters are fed to FeynHiggs up 
to the END marker. The identifier END can be cho- 
sen freely but must match exactly the one at the 
bottom of the parameter list. 

7.3. SUSY Les Houches Accord Format 

The FeynHiggs executable can also process 
files in SUSY Les Houches Accord 2 (SLHA2) 
format ^. It uses the SLHA Library [3D]. Pro- 
cessing of SLHA2 files can also be done in Library 
Mode with the subroutine FHSetSLHA. 

FeynHiggs in fact tries to read each file in 
SLHA format first and falls back to its native for- 
mat if that fails. 

7.4. Web Mode 

The FeynHiggs User Control Center (FHUCC) 
is online at feynhiggs.de/fhucc. It is a Web in- 
terface for the command-line frontend. The user 
gets the results together with the input file for 
the command-line frontend. 



7.5. Mathematica Mode 

A more powerful interactive environment is 
provided by the Mathematica interface of Feyn- 
Higgs. The MathLink executable MFeynHiggs 
must first be loaded with 

Install ["MFeynHiggs"] 

and makes all FeynHiggs routines available as 
Mathematica functions. In combination with the 
arsenal of standard Mathematica functions such 
as ContourPlot and Manipulate, even sophisti- 
cated analyses can be carried out easily. 

8. Summary: Main New Features 

Version 2.7 of FeynHiggs includes the following 
new features: 

• Inclusion of the full cMSSM two-loop 0{atas) 
corrections in highly optimized form. Extended 
interpolation options for rMSSM parts. 

• Inclusion of the full one-loop corrections to the 
hi -^ fjk decays. 

• Inclusion of NMFV corrections in the Higgs 
self-energies, decays, and low-energy con- 
straints. 

• Improvements in the A(, resummation. 

• Total Higgs production cross-sections in effec- 
tive coupling approximation, with particular 
improvements in the ggh channel. 

REFERENCES 

1. R. Peccei, H. Quinn, Phys. Rev. Lett. 38 
(1977) 1440; Phys. Rev. D16 (1977) 1791. 

2. S. Dimopoulos, S. Thomas, Nucl. Phys. 
B465 (1996) 23 hep-ph/9510220|. 

3. A. Pilaftsis, Phys. Rev. D58 (1998) 096010 
>ep-ph/9 803297 ; Phys. Lett. B435 (1998) 
88 |hep-ph /9805373 ; 

D. Demir, Phys. Rev. D60 (1999) 055006 

[hep-p h/9901389] ; 

S. Choi, M. Drees, J. Lee, P hys. Lett. B481 

(2000) 57 fhep-ph/0002287'; 

A. Pilaft sis, C. Wagner, N ucl. Phys. B553 

(1999) 3 [hep-ph/990237ll ; 



T. Hahn et al. 



7. 



10 



M. Carena, J. Ellis, A. Pilaftsis, C. Wagner, 
Nucl. Phys. B586 (2000) 92 



|hep-ph/0003180 



T. Ibrahim and P. Nath, Phys. Rev. D63 
(2001) 035009 liep-ph/0008237; Phy s. Rev. 



D66 (2002) 015005 |hep-ph/0204092| : 

S. Heinemeyer, Eur. Phys. J. C22 (2001) 

521 |hep-pli/0108059 . 

S. Heinemeyer, W. Hollik, H. Rzehak, 

G. Weiglein, Phys. Lett. B652 (2007) 300 

[arXiv:0705.0746, 

A. Brignole, G. Degrassi, P. Slavich, 

F. Zwirner, Nucl. Phys. B631 (2002) 195 

[hep-pIi/011 2177 . 

A. Brignole, G. Degrassi, P. Slavich, 

F. Zwirner, Nucl. Phys. B643 (2002) 79 

[liep-ph/0206101. 

A. Dedes, G. Degrassi, P. Slavich, Nucl. 

Phys. B672 (2003) 144 'hep-ph/0305127'. 

R. Hempfling, Phys. Rev. D 49 (1994) 6168; 

L. Hall, R. Rattazzi, U. Sarid, Phys. Rev. D 

50 (1994) 7048 hep-ph/9306309 ; 

M. Carena, M. Olechowski, S. Pokorski, 

C. Wagner, Nucl. Phys. B 426 (1994) 269 

fhep-ph/9402253 . 

M. Carena, D. Garcia, U. Nierste, 

C. Wagner, Nucl. Phys. B 577 (2000) 577 

fhep-ph/9912516;. 

S. Heinemeyer, W. Hollik, F. Merz, 

S. Pefiaranda, Eur. Phys. J. C37 (2004) 481 

fhep-ph/0403228; . 

11. S. Heinemeyer, W. Hollik, G. Weiglein, 
Comput. Phys. Comm. 124 (2000) 76 
fhep-ph/9812320. 

12. S. Heinemeyer, W. Hollik, G. Weiglein, Eur. 
Phys. J. C9 (1999) 343 [hep^ h/9812472]. 

13. G. Degrassi, S. Heinemeyer, W. Hollik, 

P. Slavich, G. Weiglein, Eur. Phys. J. C28 



(2003) 133 hep-ph/0212020 



14. M. Frank, T. Hahn, S. Heinemeyer, 

W. HoUik, H. Rze hak, G . Weiglein, JHEP 
0702 (2007) 047 jhep-ph/0611326 . 

15. T. Hahn, S. Heinemeyer, W. HoUik, 

H. Rzehak, G. Weiglein, arXiv:0709ll9Q3 

16. S. Heinemeyer, W. Hollik, G. Weiglein, Eur. 
Phys. J. C16 (2000) 139 |hep-ph/0003022 . 

17. L. Hofer, U. Nierste and D. Scherer, JHEP 
0910 (2009) 081 |arXiv:0907.5408J . 



18. A.J. Buras, P.H. Chankowski, J. Rosiek and 
L. Slawianowska, JVucL Phys. B659 (2003) 3 
|arXiv:hep-ph/021014^. 



19. http://theory.fi.infn.it/grazzini/ 
hcalculators.html 
D. de Florian, M. Grazzini, Phys. Lett. 



B 



674 (2009) 291 ;arXiv:0901.2427 , 

20. C. Anastasiou, R. Boughezal, F. Petriello, 
JHEP 0904, 003 (2009) |arXiv:081 1.3458) . 

21. R. Bonciani, G. Degrassi, A. Vicini, JHEP 
0711 (2007) 095 arXiv:0709.4227 , 

U. Aglietti, R. Bonciani, G. Degrassi and 

A. Vicini, JHEP 0701 (2007) 021 
|arXiv:hep -ph/0611266 . 

22. A. Dedes, P. Slavich, JVucJ. P hys. B 657 
(2003) 333 "hep-ph/0212132\ A. Dedes, 

G. Degrassi, P. Slavich, iVuci . Phys. B 672 
(2003) 144 he p-ph/0305127| . 

23. T. Hahn, S. Heinemeyer, F. Maltoni, 
G. Weiglein, S. Willenbrock, 
|hep-ph/0607308. 

24. K.E. Wilhams, G. Weiglein, 
arXiv:0710.5331, K.E. Williams, 

G. Weiglein, Phys. Lett. B660 (2008) 217 
[arXiv:0710.5320^. 

25. E. Braaten, J.P. Leveille, Phys. Rev. D 22 
(1980) 715. 

26. S. Heinemeyer, W. Hollik, G. Weiglein, Phys 
Rept. 425 (2006) 265 hep-ph/04122l4 . 

27. R.V. Harlander, W.B. Kilgore, Phys. Rev. 
68 (2003) 013001 |hep-ph/0304035j. 

28. T. Hahn, W. HoUik, J.I. Illana, 
S. Pefiaranda, [hep-ph/0512315, 

29. P. Skands et al^JHEP 0407 (2004) 036 
|hep-ph/03 11123 ; 

B. AUana ch et al., arXiv:080L0045| 

30. T. Hahn, hep-ph/0408283; Comput. Phys. 
Commun. 180 (2009) 1681 



D 



arXiv:hep-ph/0605049]. 



